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ABSTRACT 

In the context of cold dark matter (CDM) cosmological models, we have simulated im- 
ages of the brightness temperature fluctuations in the cosmic microwave background 
(CMB) sky owing to the Sunyaev Zel'dovich (S-Z) effect in a cosmological distribution 
of clusters. We compare the image statistics with recent ATCA limits on arcmin-scale 
CMB anisotropy. The S-Z effect produces a generically non-Gaussian field and we 
compute the variance in the simulated temperature-anisotropy images, after convolu- 
tion with the ATCA beam pattern, for different cosmological models. All the models 
are normalised to the 4-year COBE data. We find an increase in the simulated-sky 
temperature variance with increase in the cosmological density parameter f2g. A com- 
parison with the upper limits on the sky variance set by the ATCA appears to rule 
out our closed-universe model: low-f2o open-universe models are preferred. The result 
is independent of any present day observations of as ■ 
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; 1 INTRODUCTION 



A key problem in modern cosmology is the determination 
of the geometry of the universe. Classical methods of de- 
riving the cosmological parameters involve measurements of 
the redshift dependence of apparent luminosities of 'stan- 
dard candles', or the angular sizes of 'standard rulers', or 
the number densities of non-evolving objects. These meth- 
ods are expected to probe the background geometry inde- 
pendent of any structure formation. Many of these classical 
approaches have been limited by the difficulties in identi- 
fying objects, distributed over cosmic timescales, that are 
untouched by astrophysical evolution; however, it may be 
noted that progress has recently been reported using super- 
novae Type fa as standard candles (Perlmutter et al. f 998) . 

The formation of structure is dependent on the back- 
ground cosmology. Attempts have been made to exploit 
this coupling by examining the parameter space allowed 
for the cosmological constants by favoured models of struc- 
ture formation. Motivated by the discovery of anisotropy in 
the cosmic microwave background (CMB), progress in the 
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understanding of physical mechanisms responsible for the 
anisotropy spectrum and the influence of the background 
cosmology in the generation of these anisotropies have led 
to attempts at deriving constraints on the cosmological pa- 
rameters from the shape of the CMB anisotropy spectrum 
(Bond et al. 1994). Anisotropies in the CMB are usually de- 
scribed by its power spectrum: Ci = (| C; m | 2 ) represents 
the anisotropy power at multipole order I, where C; m are 
the coefficients of the spherical harmonic decomposition of 
the fractional temperature fluctuations. The density fluctu- 
ations that are hypothesized to have been generated in the 
very early universe are believed to grow via gravitational 
instabilities and give rise to the large scale structures we see 
in the present day universe. The CMB temperature fluctu- 
ations seen today on large angular scales exceeding about 
a degree are believed to be a direct consequence of mat- 
ter inhomogeneities on scales exceeding « 100 Mpc at the 
recombination epoch. The gravitational and astrophysical 
evolution, in the post-recombination universe, which led to 
the formation of galaxies and their clustering, may have al- 
tered these primary radiation anisotropies and may have 
given rise to the dominant CMB fluctuations on small an- 
gular scales. Since the generation of the CMB anisotropies 
and their appearance on the sky is intimately finked with 
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the background cosmology, the study of CMB anisotropies 
on different scales give constraints on theories of structure 
formation and the also on the parameters of cosmological 
models. 

Anisotropies in the CMB are divided broadly into two 
classes: primary and secondary anisotropies. The primary 
anisotropies are those generated at the last scattering sur- 
face. The angular power spectrum of these anisotropies, as 
seen by observers at the present epoch, depends on the pa- 
rameters of the background cosmology — the Hubble con- 
stant H , the mean matter density Q a and the constituents 
of the matter (i.e., the amount of baryonic matter, cold 
dark matter, hot dark matter etc.) and any cosmological 
constant A — and on the magnitude and spectrum of pri- 
mordial perturbations in the matter. Several authors have 
derived the marginal and joint dependencies of the cosmo- 
logical parameters on characteristic features in the CMB 
anisotropy spectrum. For example, at large and intermedi- 
ate angular scales (I < 500), for Qo < 1, the Doppler peaks 
in the CMB anisotropy spectrum are shifted towards larger 
/ values (see, for example, Kamionkowski & Spergel 1994 or 
Lineweaver & Barbosa 1998). The ionization history modi- 
fies the radiation spectrum and serves to attenuate power at 
high-order multipoles. Primary anisotropies at small angu- 
lar scales (I > 500) are expected to be relatively damped in 
most models of structure formation due to the thickness of 
the last scattering surface and the diffusion damping of sub- 
horizon scale baryon fluctuations in the pre-recombination 
era; anisotropies at even larger I may be critically dependent 
on the reionisation history. 

The secondary anisotropies arise due to interactions be- 
tween the growing matter perturbations and the CMB pho- 
tons as they travel from the last-scattering-surface to the 
present time. Arcmin-scale anisotropies may be generated 
due to second-order mode coupling between density pertur- 
bations and bulk velocities (Hu, Scott & Silk 1994; Persi 
et al. 1995). Separately, decrements in the CMB sky may 
be generated owing to the inverse-Compton scattering of 
background photons as they pass through concentrations 
of hot gas associated with clusters of galaxies along the 
line of sight; this phenomenon is referred to as the Sun- 
yaev Zel'dovich (S-Z) effect. A dominant contribution to the 
anisotropy power at large / may come from quasar-ionised 
hot gas bubbles (Aghanim et al. 1996). 

Among the various effects that may give rise to CMB 
anisotropy on small angular scales, we focus on the S-Z 
effect in this work. In section 2, we discuss some details 
of the S-Z effect and in section 3, we discuss some of the 
motivations for attempting to constrain cosmological mod- 
els using the S-Z effect. In section 4, we present the for- 
malism we have adopted for deriving the expectations for 
CMB anisotropy on small angular scales for specific mod- 
els of structure formation. The results of simulations of sky 
temperature anisotropy are compared with ATCA limits on 
arcmin-scale anisotropy to derive constraints on the cosmo- 
logical density parameter Q Q . 



2 THE SUNYAEV-ZEL'DOVICH EFFECT 

The hot intracluster gas inverse-Compton scatters the CMB 
photons propagating through the cluster medium, and the 



energy transfer in this interaction between the hot electrons 
and CMB photons results in a distortion to the CMB spec- 
trum (Sunyaev & Zel'dovich 1972). The integral of the elec- 
tron pressure along any line-of-sight through the cluster de- 
termines the magnitude of the distortion in that direction, 
this is quantified in terms of the Compton y-parameter: 



y = /« 



h(T e T<y 

2 
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where k is the Boltzmann's constant, T e and T 7 are, re- 
spectively, the kinetic temperature of the electron and the 
thermodynamic temperature of the incident black-body ra- 
diation, m e is the electron rest mass, c is the velocity of light, 
n e is the electron number density and ot is the Thomson 
scattering cross-section. 

The S-Z effect manifests itself as a change in sky bright- 



<5v = yj v (x), 



(2) 



(3) 



towards the cluster with respect to the mean CMB intensity, 
x is a dimensionless frequency parameter defined to be 

hv 

where h is the Planck's constant, v is the observing fre- 
quency and T is the CMB temperature at the present 
epoch: T — 2.73 K. The frequency dependence of the change 
in sky brightness, owing to the S-Z effect, is given by 



2(fcT 



tanh (x/2) 



(4) 



The spectral function has a unique shape owing to the fact 
that although there is a net transfer of energy from the hot 
electron gas to the radiation, inverse-Compton scattering 
conserves the total number of photons and the mean photon 
energy is increased. The result is that there is a net diffusion 
of photons upwards in energy. Assuming the cluster electron 
temperature to be about 10 keV, the net distortion observed 
is expected to be zero at about 222 GHz (1.35 mm wave- 
length), the distortion is expected to result in a decrement 
in the brightness temperature towards the cluster at lower 
frequencies and the cluster will be seen as a positive source 
at higher frequencies. 

At wavelengths longward of 1.35 mm, a cluster appears 
as a negative source on the sky with respect to the mean 
CMB background intensity: the flux density S v due to the 
integrated S-Z effect over the sky area of a cluster is 



S v {x) = 



Dl(z) 



1 



dV 5-n e crT. 

m e c 2 



(5) 



We have assumed that T e ^> T 7 . The integral extends over 
the cluster volume and D a (z) is the angular-size distance to 
the cluster. It is clear from the equation above that the net S- 
Z effect due to a cluster is dependent only on the total mass 
of hot gas in the cluster and a density-weighted mean tem- 
perature. The S-Z flux density from a cluster diminishes as 
the square of its angular-size distance: because the angular- 
size distance of objects at cosmological distances, z <; 1, 
saturates to a limiting value or even decreases with increas- 
ing redshift, the S-Z flux from a cluster does not rapidly di- 
minish with increasing cosmological distance. However, the 
distributions of cluster gas temperature and electron density 



© 0000 RAS, MNRAS 000, 000-000 



Structure Formation and SZE 3 



may be functions of redshift; therefore, the expected S-Z flux 
density from a cluster may be redshift dependent. 

It may be noted here that y is a dimensionless parame- 
ter and S„ has units of erg s" 1 cm 2 Hz -1 : this flux density 
is usually expressed in Jansky (Jy), defined as 1 Jy = 10 -23 
erg s _1 cm 2 Hz -1 . 



3 MOTIVATION FOR STUDYING THE S-Z 
EFFECT 

The Sunyaev-Zel'dovich effect provides a rather nice and 
complementary approach to the traditional methods - 
which use X-ray temperature and X-ray luminosity evolu- 
tion in clusters of galaxies — of studying the evolution of 
the mass function of collapsed objects (Bartlett 1997). The 
evolution in the abundance of clusters of galaxies is sensitive 
to the mean matter density fi and, consequently, is a useful 
constraint on cosmological models. The traditional method 
derives its faith from numerical simulations which show that 
there is a tight correlation between the virial mass and the 
emission- weighted X-ray temperature in clusters of galaxies 
(Evrard 1990, 1996). However, the traditional method has 
a major disadvantage: because of 'cosmological dimming', 
the surface brightness of distant X-ray sources falls off as 
(1 + z)~ 4 , and for this reason, obtaining samples of clusters 
at cosmological distances can be observationally challeng- 
ing. In the case of the S-Z effect in cosmological clusters, the 
decrement in brightness in lines of sight towards clusters of 
galaxies has the distinct advantage of being independent of 
the distance to the cluster. The S-Z flux density from a clus- 
ter will diminish with distance to the cluster as the square 
of the angular-size distance; this is in contrast to X-ray flux 
densities from clusters which diminish as the square of the 
luminosity distance to the cluster. Moreover, the integrated 
S-Z effect due to any cluster, i.e., the flux density decre- 
ment, is proportional to the total hot-gas mass times the 
particle-weighted temperature. Consequently, the detection 
of any cluster is independent of the gas' spatial distribution 
(assuming the cluster is unresolved). If the observations re- 
solve clusters, particularly at lower redshifts, the observed 
sky temperature distribution will be sensitive to the temper- 
ature structure within the clusters; once again this may be 
contrasted with X-ray emission images of cluster gas distri- 
butions which are sensitive to the gas density distribution. 

The S-Z effect may, therefore, be used as a tracer of 
clusters and other massive hot gaseous objects at cosmolog- 
ical distances. Because the redshift evolution of clusters is a 
sensitive probe of cosmology theory, the number counts of 
S-Z sources may be used to infer their cosmological abun- 
dances and thereby deduce parameters relating to the back- 
ground cosmological model and the structure formation the- 
ory (see, for example, Blanchard & Bartlett 1998; Oukbir & 
Blanchard 1992, 1997). Motivated by this reasoning, several 
authors have, in recent years, made calculations of the ex- 
pected S-Z source counts and their redshift distributions. 
These have been related to the cosmological parameters as 
well as the evolution in the intra-cluster medium (Bartlett 
& Silk 1994; Barbosa et al. 1996; Colafransesco et al. 1997). 

As discussed in section 2, clusters of galaxies may be 
'visible', owing to the S-Z effect, as sources in the sky with 
a negative flux density at wavelengths longward of 1.35 mm 



and with a positive flux density at shorter wavelengths. This 
S-Z effect has been imaged to-date towards several clus- 
ters (see Birkinshaw 1999 for a recent review). We may 
consider distant cosmological clusters — observed at wave- 
lengths longward of 1.35 mm — as a sky distribution of 
sources with negative flux density. Observational sensitivity 
in modern radio telescopes is reaching values close to that 
required for detecting the S-Z effect from cosmological clus- 
ters towards 'blank fields' where no obvious clusters are seen 
in either their optical or X-ray emission. There have been 
claims in the literature of the detection of radio decrements 
(thought to be due to the S-Z effect) in sensitive images 
made of 'blank' sky fields (Jones et al. 1996; Richards et al. 
1996); however, sensitive observations of several other fields 
with arcmin resolution have failed to detect any decrements 
or CMB anisotropies (Subrahmanyan et al. 1998, 1999). 

It may be noted here that because the S-Z effect has a 
generic non-Gaussian temperature distribution, it could be 
detected in a sky image by estimators sensitive to skewed 
variance. Many other secondary contributors to the temper- 
ature anisotropy have a Gaussian distribution in amplitudes; 
therefore, the negative-skewed-nature of the S-Z effect may 
be useful in distinguishing it from other sources, foregrounds 
and instrument noise. 



4 SIMULATIONS 

4.1 The cosmological distribution of clusters of 
galaxies 

As discussed earlier, large-angle anisotropies in the mi- 
crowave background radiation may be traced back to its 
generation from small-amplitude primordial density fluctu- 
ations in the early universe. The matter perturbations have 
since grown due to gravitational instabilities into the large 
scale structure we see at the present epoch. In addition, as 
the primeval CMB radiation propagated to us through grow- 
ing matter inhomogeneities, astrophysical couplings give rise 
to secondary anisotropies on predominantly small angular 
scales. The power spectrum P(k, z eq ) of the fractional den- 
sity fluctuations, at the redshift z — z eq when the energy 
densities in matter equalled that in radiation, may be re- 
lated to the primordial power spectrum P p (k) by: 

P(k,z eq )=T 2 (k)x P p (k). (6) 

The matter transfer function T(k) describes the process- 
ing of the initial density perturbations during the radiation 
dominated era (Padmanabhan 1993): the modifications to 
the shape of the perturbations would depend on the nature 
of the perturbations and the candidate dark matter. There- 
after, during the matter-dominated era, the dominant dark 
matter (DM) perturbations experience equal growth on all 
scales and P(k, z) grows retaining its shape. The growth 
rate, described below by the growth function, varies with 
cosmic time and depends on the mean matter density . 
Models of structure formation are usually characterised by 
specific P(k,z eq ) shapes and normalizations. 

For the purpose of our simulations, we wish to relate 
the number density of collapsed objects of different masses, 
at different cosmic epochs, to the initial density contrast on 
different scales, i.e., the initial matter spectrum P(k,z eq ). 
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It is assumed, in most theories of structure formation, that 
the initial small-amplitude density fluctuations are Gaussian 
random, i.e., the amplitudes are Gaussian distributed and 
that the fluctuations on different modes are uncorrelated 
(random phase fluctuations). Structure is believed to form 
from the growing perturbations hierarchically, with smaller- 
scale fluctuations collapsing first and larger scales later. The 
mass and redshift distribution of the number density n(M, z) 
of collapsed objects may then be computed using the Press- 
Schechter formalism (Press & Schechter 1974; see Padman- 
abhan & Subramanian 1992 for a tutorial): 



n(M,z)dM = «/-Mj/ c 

y 7r m 

where 



d\no(M,z) 



dlnM 



S c D g (0) 



-v 2 J2 



dM 

M ' 



(7) 



(8) 



Vc{M,Z) ~ <r(M,z) - ao(M)Dg(z) 

In these equations, (p) denotes the mean co-moving matter 
density in the universe: (p) = 30. o H 2 /(8ivG), where is the 
density parameter, H a — lOO/i km s~ x Mpc" 1 is the Hub- 
ble constant and G is the gravitational constant. <j(M, z) 
denotes the rms fluctuations, at redshift z, in the fractional 
density contrast in the matter when smoothed to a mass 
scale M = ^TrR 3 {p), where R is the comoving radial length 
scale of the smoothing function. a (M) is the rms density 
contrast at the present epoch. 

The density contrast in an overdensity, computed us- 
ing linear theory, at the epoch when the collapsed object 
is deemed to have formed, is denoted by S c . We adopt the 
collapse of a spherical top-hat overdensity (Peebles 1980) as 
a valid model for the description of the dynamical evolu- 
tion of the peaks in the Gaussian density distribution; the 
collapsed object is deemed to have 'formed' at the cosmic 
time t c , which is approximately twice the cosmic time t m 
at which the overdense region attains maximum expansion 
radius. In a universe with Q = 1.0, S c = 1.68. It has been 
shown that S c varies by at most ~ 4 per cent for a range of 
O from 0.1 to 1 (Lacey & Cole 1993). Therefore, we have 
chosen to adopt a constant value 1.68 for 8 C . 

fn this Press-Schechter formalism, collapsed objects are 
identified on the basis of their overdensity assuming linear 
growth of perturbations. D g , in equation 8, quantifies the 
growth factor of the density perturbations from the epoch 
of matter- radiation density equivalence (z = z eq ) to any 
later epoch z in the matter dominated era. In the absence 
of free streaming, the growth function is given by (Peebles 
1980, Heath 1977) 



D g (z) = ^(l + z eq )g(z) 
where 



1+z' 



dz , 



g 2 (z) = fi (l + zf + (1 - fi„ - 00 (1 + zf + Ov. 



(9) 



(10) 



Closed form expressions are available in Weinberg (1972), 
Groth & Peebles (1975) and Edwards & Heath (1976) for 
universes with and without a cosmological constant A; we 
have used these expressions for the purpose of our simula- 
tions. 

The rms amplitude of the mass fluctuations at any red- 
shift z: a(M,z), when smoothed with a spherically sym- 
metric window function of characteristic co-moving radius 



7?, may be computed from the the matter power spectrum 
P(k, z) using the relation: 



a 2 (M,z) 



t^P(K*)\w R (k)\\ 



(ii) 



where Wn{k) is the Fourier transform of the corresponding 
real space window function and, as before, M — ^ivR 3 {p}. A 
spherical top-hat form with radius R is usually adopted for 
the window function and this corresponds to a Fourier-space 
window function: 

W R (k) = -^3 (sm(kR) - kRcos(kR)). (12) 
kR 

The Press-Schechter formalism has been verified by 
Tozzi & Governato (1998) against the results from numeri- 
cal n-body dynamical simulations, particularly for the mass 
range encompassing clusters of galaxies. 

4.2 P(k) and its normalization 

We have adopted a power-law form for the initial primordial 
matter spectrum with P p (k) ~ k n , where n is the index of 
the primordial power spectrum; this index n equals unity 
for Harrison-Zel'dovich scale-invariant spectra. The power 
spectrum at any redshift z in the matter dominated era may 
then be written in the form 



k 3 . (ck\ 3+n 



S 2 H T 2 (k)D 2 g (z)/D 2 g (0), 



\H J r '" s "'" sr " (13) 

where the D g 's are the growth factors defined in section 4.1. 
8h represents the amplitude normalization of P(k, z) and 
is defined as the amplitude of perturbations on the horizon 
scale at the present epoch z = 0. 

The large-angle CMB anisotropics discovered by 
COBE (Smoot et al. 1992) are believed to have been gen- 
erated by processes at the last scattering surface. Precisely 
for the reason that the matter fluctuations were of small 
fractional amplitudes at the epoch t they generated these 
anisotropics, and because the scale length of the pertur- 
bation modes well exceeded ct, the coupling physics be- 
tween the matter perturbations and radiation was linear and 
did not involve any astrophysical interactions. The CMB 
anisotropies detected by COBE may, therefore, be related 
to the primordial matter spectrum fairly exactly without all 
the uncertainties associated with the astrophysics involved 
in the formation of collapsed baryonic objects at late epochs. 
For this reason, we have chosen to normalize the P(k) with 
the COBE results. 

We have used the fitting formulae of Bunn & White 
(1997) to normalize the power spectrum to the COBE- 
DMR measurements. For an open universe with vanishing 
cosmological constant, with no CMB anisotropies coming 
from gravitational waves, we use the normalization 



c i n r .. 1n -5 o -0.35-0.191nn -0-17ri -n-0.14fi 2 

Oh = 1-95 x 10 i! u e 



(14) 



n — n — 1 and the normalization is valid for 0.7 < n < 1.2. 
The fit is also valid over the range 0.2-1.0 in O- 

As mentioned earlier, the shape of P(k, z) in the mat- 
ter dominated era is completely determined by the trans- 
fer function T(k) and the assumed form for the primordial 
spectrum. With the choice of a power-law form for P p (k), 
the transfer function is critical in relating the power on small 



© 0000 RAS, MNRAS 000, 000-000 



Structure Formation and SZE 5 




\ 



- I i L i I I ^- I , i I i l_ 

1 2 4 1 2 4 
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Figure 1. Abundance of objects with mass M = 10 12 M (continuous line), M = 10 13 M (dashed line), M = 10 14 Mq (dot-dashed 
line) and M = 10 15 Mq (dotted line). A CDM cosmology with h = 0.65, A = and Q B = 0.019/i~ 2 is adopted; the four plots arc 
labelled with the assumed values of Q . The abundance N(M, z) is expressed as the number of objects per Mpc 3 and logio(M X N(M, z)) 
is plotted versus redshift z. 



wavenumber modes — which is fixed by the adopted COBE 
normalization — to those on large k modes at which we have 
clusters of galaxies. There are several fitting forms avail- 
able in the literature for T(k) in the context of the cold 
dark matter (CDM) dominated universes. However, different 
parameterizations of transfer functions can differ by large 
amounts (Peacock & Dodds 1993) and, moreover, to obtain 
better accuracy the effect of baryon damping must be in- 
cluded (Hu & Sugiyama 1996; Ma 1996; Eisenstein & Hu 
1999). We, however, work within the context of low-baryon- 
density universes and, therefore, include the effect of non- 
zero baryon content by adopting a modified 'shape factor' T. 
In this work, we have adopted the fitting function provided 
by Bardeen et al. (1986): 

In (1 + 2.34g) 
-LCDM(q) = z-z-. x 



2.34g 



[l + 3.89q + (16.1g) 2 + (5.46g) 3 + (6.71g) 4 ] 1/4 (15) 



where q = k/hT and F — Qohexp (— Q,b —Qb/^o)- is 
the baryon density parameter and is the ratio of the mean 
density of baryons in the universe to the critical density 
Pc = 3H 2 /(8ttG). 

The Press-Schechter formalism described in section 
4.1, along with the COB .©-normalized P(k) defined above, 
have been used to calculate the number densities of col- 
lapsed objects. We choose to identify the condensates as 
groups/clusters of galaxies and have computed their abun- 
dances as a function of cluster mass and redshift for a set of 
CDM cosmologies differing in their Q . We have not consid- 
ered the effect of including a cosmological constant in this 
work and have put A = 0. The computations of cluster distri- 
butions n(M, z) have been carried out over logarithmic bins 
in the mass range 1O 13 -1O 16 M . The distribution in mass is 
computed at redshifts spaced at intervals of Az = 0.1. Plots 
of n(M, z) distributions are shown in Fig. 1. 

The theory predicts that the cosmic abundance of mas- 
sive collapsed objects are extremely sensitive to the ampli- 
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tude and slope of the primordial power spectrum and also 
on the growth factor. In a flat cosmology with fl = 1, 
P(k) would grow as (1 + z)~ 2 ; in an open universe the 
growth approximately follows this evolution down to red- 
shift z w (flJT 1 — 1) and is stunted thereafter. Therefore, 
in these models, which are all COBE normalized, the to- 
tal growth up to the present time will be greater in models 
with larger £7 . Consistent with this reasoning, it may be 
seen in Fig. 1 that the larger Q a models are dynamically 
more evolved: they have larger numbers of 10 15 Mq objects 
at z — and in these models the abundances of lower mass 
objects are declining with cosmic time as they are incorpo- 
rated into larger mass objects. 

It may be seen from Fig. 1 that the decline, with red- 
shift, in the abundance of high-mass (10 15 M ) collapsed 
objects relative to the abundance at z = is steepest in the 
case of models with lower Q a . This is the opposite of the ex- 
pectations for models which are normalized at the present 
time to, for example, a specific value of as, which is the rms 
mass fluctuations when the present day universe is smoothed 
using a window function with a radius R — 8/iT 1 Mpc. Al- 
though the growth function evolves slower at low redshifts in 
models with low £l , the abundance of massive collapsed ob- 
jects evolves more strongly in the case of the low fi models 
because in these models the rms mass fluctuations a(M, z) 
is itself low at z = 0. 

The Gaussian characteristic of the Press- Schechter mass 
function is evident in the distributions: at any redshift, there 
is a rapid exponential decline in the number density of ob- 
jects in larger mass bins. 



4.3 Simulated sky images 

Our simulation codes compute the abundances n(M, z) in 
redshift slices in the redshift range z — 0-5. We make 
the conservative assumption that collapsed objects in the 
restricted range 10 13 -10 16 Mq alone represent clusters of 
galaxies containing hot gas. Adopting a model, described 
below, for the spatial distribution of hot gas in the poten- 
tial wells of these clusters, we simulate sky images of the 
S-Z effect expected owing to these clusters. The sky was 
simulated by separately computing the contributions from 
different redshift slices along the line of sight. In any red- 
shift slice, the clusters are assumed to be Poisson-random 
distributed on the sky. We have simulated square patches of 
sky consisting of 256 pixels along each side: the pixels were 
chosen to be 10 arcsec square making the total image size 
42'40". The redshift slices were of size Az = 0.1. 

The mean number of clusters, A, with mass in the range 
M-(M + 5M), is related to the Press-Schechter number den- 
sity by: 

A = n(M, z) x 8M x volume corresponding to a pixel. (16) 

The comoving volume of any pixel in a slice at redshift z is 
calculated as the product of the comoving area (AZ) 2 covered 
by the square pixels of angular size A6 = 10 arcsec and 
the comoving line-of-sight distance As corresponding to the 
redshift slice Az. These are given by: 



Al 



2cA6 



As = 



dz 



H (1 + z)y/l + 0, o Z 



(17) 



HM'l(l + z) 



{D. Z + {Qa - 2) [v^o + 1 - l] } . (18) 



We have associated every collapsed massive object, with 
mass exceeding 10 13 Mq, with cluster gas. In order to model 
the spatial distribution of the intra-cluster gas, the isother- 
mal /3 model has been adopted, in which the cluster gas is 
modelled as being spherically symmetrical, centred in the 
gravitational potential of the cluster and shock heated to a 
temperature corresponding to the infall energy. The varia- 
tion in gas density with radial distance r is assumed to be 
given by 



p(r) = p [l + (r/r c ) 2 ] 



-3/3/2 



(19) 



where p a is the central density and r c is the core radius. 
The value of (3, inferred from the X-ray surface-brightness 
profiles observed in clusters of galaxies, is believed to be in 
the range 0.5-0.9 (Markevitch et al. 1998, Jones & Forman 
1984). Herein, for simplicity, we adopt a value of | for (3. 

Assuming that the dynamical collapse of the cluster is 
self similar, the scaling of the core radius r c of collapsed 
objects with mass and redshift has been derived to have the 
form (Colafrancesco 1997, Kaiser 1986): 

, 1.3/i _1 Mpc 1 
r c (ft ,M, z) = — — — r X 



M 



o (1 + z) 
A(ft = M = 0) 



-, 1/3 



10 15 h-iMe Sl A(n ,z) 



(20) 



In this parametrization, the non-linear density contrast at 
the epoch of virialisation is A = p/p where p is the density 
contrast and p is the mean background density. For the case 
of an open universe, 



18tt^ 



1 



(21) 



and 



where t v and z v are the cosmic time and redshift corre- 
sponding to the epoch of virialisation. In the case of a flat 
universe, A ^ 178 (Peebles 1980; Colafransesco et al. 1997). 
In the above parametrization, p is an adjustable free pa- 
rameter that may be selected to fit observations of nearby 
clusters. If r m denotes the maximum radial extent of the 
cluster, p — (r m /r c ) and is a measure of the concentration 
of the cluster mass: small values of p model the cluster gas as 
centrally concentrated where as larger values of p spread the 
gas mass away from the core. Values of p « 6 appear to bet- 
ter model the parameters of nearby clusters; a comparison 
of the model with the observed parameters for the Coma 
cluster is given below. We propose to use the simulations 
for the purpose of predicting the sky variance as observed 
by the ATCA; these observations are made with arcmin res- 
olution and will be less sensitive to S-Z effects from clusters 
whose gas distribution is extended. The ATCA observations 
couple somewhat better to models with smaller p; therefore, 
we have adopted the conservative choice of p = 8. It may 
be noted that sensitive X-ray imaging of the gas in nearby 
clusters have shown gas extending out to at least eight core 
radii. 

Although recent observations indicate that intracluster 
gas has temperature structure (Markevitch et al. 1998), we 
have modelled the clusters as being isothermal. It follows 
from the assumption that clusters form from self-similar col- 
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lapse that the temperature T « and a good approxima- 
tion for the temperature of the intra-cluster gas in clusters 
is (Colafrancesco et al 1997) 

2/3 

6.7 x W r (l + z) 



M 



10 w h~ 1 M, 



o 



r2 A(n ,2) 



-, 1/3 



A(f2 = l,z = 0) 



K. 



(22) 



The relation between the mass and the temperature in equa- 
tion(22) agrees well with the recent M-T parametrisation as 
suggested by Bryan & Norman (1998). We have assumed 
that the hot gas detected in its X-ray emission is respon- 
sible for any S-Z effect as well (Colafrancesco et al. 1994; 
Blanchard et al. 1992; Vittorio et al. 1997; Kaiser 1986). 

Only the baryonic matter in collapsed clusters gives the 
S-Z effect. White et al. (1993) have estimated the baryon 
content in clusters to be 
M b 



> 0.009 + 0.050h~ 3/2 , 



(23) 



where the first contribution is due to the galaxies and the 
second is due to the intra-cluster gas. Primordial nucleosyn- 
thesis predicts a universal baryon abundance Q, B ~ 0.019/i~ 
(Buries et al. 1999). Adopting the Hubble parameter h = 
0.65, it is seen that the baryon gas-mass fraction in nearby 
clusters is about 0.1 and is a factor of two greater than the 
universal baryon abundance. The Einstein Medium Sen- 
sitivity Survey (EMSS) data appears to indicate a decline 
in the abundance of X-ray luminous clusters with redshift 
(Henry et al. 1992); this may be parametrized as a hot-gas 
fraction that evolves as t 1A with cosmic epoch. It may be 
noted that adopting a gas fraction that declines with red- 
shift is a conservative assumption because it reduces the 
predicted S-Z effect. We have adopted a parametrization 

-i 0.2 



M 



0.050fr~ 3/2 



M u 



lO 15 ^- 1 M 



© 



(24) 



with q — 1.4. Mgas represents the hot gas mass in the clus- 
ter and Mtot represents the total mass in the collapsed ob- 
ject. The exponent 0.2 in equation(24) is from a fit in Co- 
lafrancesco et al (1997). We have assumed that the intra- 
cluster gas has a primordial composition, with helium and 
hydrogen atoms in the number ratio 1:10, and is fully ion- 
ized. 

The cluster mass distribution, as also the intra-cluster 
gas, are modelled as truncated at a maximum radius r m = 
p x r c . We define the 'impact parameter' b as the projected 
distance between the cluster centre and any line of sight 
through the cluster. It follows that the S-Z decrement in the 
Rayleigh- Jeans portion of the spectrum, in units of temper- 
ature, expected at any sky position b Si r ma x is given by 



AT : 



-TcMB 



2kT 

m e c 2 



n e <jT<ll, 



(25) 



where Tcmb is the average brightness temperature of the 
CMB. The integral is computed along the line of sight 
through the cluster at the impact parameter b. 

The assumption that the cluster gas is of primordial 
composition implies that the total number of hot electrons 
in the cluster is iV e = (18/14)(M gas /a.m.u). The S-Z decre- 
ment may be written in the form 



AT/Tcmb 



2kT 

m e c 2 



-a T x 



tan -1 ( 



\J~r 



-b 2 



)N e r c 



4tt ^ 



n + b 2 [r c r„ 



. r 2 tan -l(I^)] 



(26) 



The Coma cluster at a redshift 0.0235, with a total mass 
about 2 x 10 15 Mq , is observed to have X-ray gas with a 
temperature 8-9 keV and a core radial size of 10.5 arcmin 
(Silverberg et al. 1997, Herbig et al. 1995). The Compton- 
y parameter has been measured to be 9 x 10~ 5 towards its 
centre (Herbig et al. 1995). In our simulations, for the choice 
p — 6 and adopting a value of 0.096 as the baryon gas frac- 
tion, we find that a collapsed object of mass 2 x 10 15 Mq 
at redshift z — 0.0235 — corresponding to the Coma cluster 
parameters — yields values of 10.3-11.0 arcmin for the core 
size, 8-9 x 10 7 K for the gas temperature and the central S-Z 
Compton-j/ decrement is in the range 7-8.5xl0 -5 . 

We have accumulated the S-Z effect from clusters in red- 
shift slices up to a maximum redshift of 5. In each redshift 
slice, the pixels were populated by collapsed point masses 
in a Poisson random fashion; the expectation that any pixel 
was populated by objects in any mass bin was governed by 
equation 16. We then substituted our cluster gas model for 
every mass point: the S-Z effect owing to each cluster is dis- 
tributed over several pixels surrounding that at which the 
point mass was located. The variance of the temperature 
fluctuations in the cumulative S-Z effect images were com- 
puted. 

We have compared the results of the simulations with 
the observational limits on arcmin-scale CMB anisotropy 
set by the ATCA experiment at 8.7 GHz (Subrahmanyan et 
al. 1993, 1998, 1999). These observations have been made 
with a Fourier synthesis array. The sky is viewed by the 
telescope as attenuated by a primary beam that has a full 
width at half maximum of about 6 arcmin and the synthetic 
images represent the convolution of the visible sky with a 
synthetic beam pattern. The ATCA experiment sets upper 
limits on any residual sky fluctuations apart from the in- 
strument noise and confusing foreground sources. Assuming 
'flat-band' CMB fluctuations, the ATCA limit corresponds 
to Qfiat ^ 25 piK and, for their telescope filter function, 
this corresponds to an upper limit of 21 ujy beam -1 on the 
flux-density fluctuations on the sky. We have convolved the 
images resulting from our simulations by a beam that is the 
product of the ATCA synthetic beam and primary beam. 
Sample images are shown in Fig. 2 corresponding to CDM 
models with f2 = 1.0 and 0.6. Images are displayed using a 
1-arcmin Gaussian (normalized to give unity volume) as the 
convolving beam and separately using the ATCA product 
beam (normalized to have a peak unity) as the convolving 
function. It is seen that the extended clusters, seen in the 
flo = 1.0 simulation, are resolved by the ATCA beam: the 
extended S-Z structures will not be detectable by the ATCA 
imaging. It may be noted that before the convolution, the 
image pixels are in units of temperature (K); following the 
convolution, the pixel intensities are in units of flux density 
(Jy beam -1 ). The predictions for the image variance have 
been made for different plausible cosmological models. 
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sCDM65 sCDM65 




arcmin arcmin 



Figure 2. Sample images. The upper pair of images corresponds to CDM cosmology with Q a = 1.0; the lower pair corresponds to 
Qo = 0.6. The cluster gas distribution has been parametrized with p = 8 and the baryon gas fraction is assumed to evolve with redshift 
as described by equation 24; the Hubble parameter h = 0.65 has been adopted. The images on the left have been convolved with a 
Gaussian beam of FWHM 1 arcmin; the images on the right have been convolved by the ATCA product beam. Contours at (—3, —2, 
— 1, 1, 2, 3, 4, 6, 8, 12, 16, 24) X — 30 /^K for the images on the left and —30 fijj beam -1 for the images on the right. 



Table 1. Statistics for different cosmological models. 



Model 


n 


h 


C8 


gas fraction^] 


sCDM 


1.0 


0.50 


1.128 


parametrized 


sCDM(high h) 


1.0 


0.65 


1.567 


parametrized 


oCDMl 


0.8 


0.65 


1.326 


parametrized 


oCDM2 


0.6 


0.65 


1.004 


parametrized 


oCDM3 


0,1 


0.65 


0.605 


parametrized 


sCDM(variantl) 


1.0 


0.50 


1.128 


no evolution 


sCDM(variant2) 


1.0 


0.50 


1.128 


M tot ~ n B 


sCDM(p=6) 


1.0 


0.50 


1.128 


parametrized 



mean y Image rms Image rms 

Xl0 — 5 before convolution after convolution 

(/iK) (A*Jy beam -1 ) 

0.61 81 28 

1.9 313 58 

0.46 60 23 

0.18 38 15 

0.05 9 3 

1.9 207 72 

2.4 245 79 

1.16 98 37 



a The parametrization of the gas fraction is described in equation 24. 
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5 RESULTS AND DISCUSSION 

Results of the simulations are given in Table 1 for a range 
of model parameters. The image variances have been listed 
both before and after convolving with the ATCA product 
beam. The standard CDM model with h — 0.5, fl a = 1.0 
and with the hot gas fraction parametrized as described in 
equation 24, is expected to result in an image rms which 
exceeds the ATCA limit. Increasing the Hubble parameter 
to a more likely value of h — 0.65 increases the expectations 
for the image rms and this 'high-/i' CDM model is rejected 
with greater confidence. Considering open universe CDM 
models with S7 < 1, the expected image variance decreases 
with decreasing Q . Low fi < 0.8 open-universe models are 
allowed by the ATCA limits. 

If we assume that q = and that the cluster gas frac- 
tion does not diminish with redshift, high-redshift clusters 
contribute to the net S-Z sky fluctuations. This is seen in the 
table where the image rms rises from 28 to 72 /uJy beam -1 
when the no-evolution assumption is made for the standard 
CDM model. We have also simulated images corresponding 
to the case where the gas fraction is a constant and equal to 
the mean baryon density predicted by nucleosynthesis. This 
gas model leads to predicted rms values which are not very 
different from that given by the q = parametrized model; 
this indicates that the lower mass (M to t ~ 10 13 ~ 14 Mq) ob- 
jects dominate the variance contribution in the no-evolution 
case. 

Simulations of images corresponding to the choice p = 

6 leads to enhancements in the expectations of the image 
rms values confirming that our earlier choice (p = 8) was 
relatively conservative. 

We have computed the mean Compton-y parameter 
from the images for the different models and these are also 
listed in Table 1. The COBE-FIRAS experiment has set 
an upper limit of y < 1.5 x 10 _s (Fixsen et al. 1996). It may 
be noted that our sCDM model with h = 0.65 is indepen- 
dently rejected by these limits. Likewise, the models with no 
evolution in gas fraction are also disallowed by the FIRAS 
results. 

The value of as has also been computed for each model 
and listed in Table 1. The standard CDM models normalized 
to COBE are known to be incompatible with estimates of 
as derived from the local X-ray luminosity function or other 
local measures of galaxy clustering (which prefer low values 
of as about 0.6). The ATCA results independently provide 
evidence in favour of a low-fi universe. It must be empha- 
sised once again that the sky fluctuations due to SZE is 
generically non-gaussian in nature and a precise elimination 
of models have to take this into account. 

An important difference between this work as compared 
to previous predictions for CMB anisotropics from S-Z ef- 
fects is that we have normalized our matter power spectrum 
to the COBE anisotropy results. All the previous predic- 
tions that we are aware of have normalized their models 
to observations of the present day clustering in galaxies or 
to X-ray luminosity functions; consequently, other workers 
have essentially normalised their matter power spectra to 
as. This makes it possible for all the models to be consis- 
tent with the observational estimates of present day cluster 
abundances. Our choice of COBE normalisation results in 
a as that varies across the models. This results in a depen- 



dance on Q which is opposite of what usually found. The 
difference may be understood in the following way : if the 
as is held constant across cosmological models, as previous 
workers have done, varying fi changes the growth func- 
tion and consequently the abundance of clusters at redsifts 
z > will decrease with increasing fl . On the other hand, 
when the matter power spectrum is normalised to COBE- 
DMR data, varying fi across models alters the shape of 
the matter spectrum and the normalisation in addition to 
the growth function. In this case, increasing 57 results in an 
increase in cluster abundances at all redshifts. Our choice of 
normalization results in a as that varies across the models. 



6 CONCLUSION 

We have considered cosmological models composed of cold 
dark matter and baryons (and no cosmological constant), 
having an initial scale-invariant spectrum of adiabatic per- 
turbations. We have used the Press-Schechter formalism to 
generate the distribution of clusters. We have normalised the 
rms mass fluctuations to COBE— DM R data. The Sunyaev- 
Zel'dovich decrement is calculated for each cluster adopting 
a model for the density profile of the cluster gas. We have 
simulated blank sky surveys for the S-Z effect owing to a 
cosmological distribution of clusters. Finally, we have pre- 
dicted the expectations for the variance in the background 
sky both before and after convolving with the ATCA beam. 
Based on a comparison with the upper limits set by the 
ATCA on CMB anisotropics on arcmin scales, we conclude 
that COB-E-normalized CDM models of the universe with a 
high density parameter (fi > 0.8) are rejected. This result 
is independent of any present epoch measures of as- 
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